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We classify entanglement singularities for various two-mode bosonic systems in terms of catastro- 
phe theory. Employing an abstract phase-space representation, we obtain exact results in limiting 
cases for the entropy in cusp, butterfly, and two-dimensional catastrophes. We furthermore use 
numerical resuhs to extract the scahng of the entropy with the non-hnearity parameter, and discuss 
the role of mixing entropies in more complex systems. 
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For a large number of quantum critical systems, criti- 
cality manifests itself as a peak, or indeed a divergence, in 
the entanglement of the ground state. Systems in which 
this behaviour has been observed include spin-1/2 ferro- 
magnetic chains in a magnètic field driven, dissipative 
large- j pseudo-spin models the Lipkin-Meshkov-Glick 
Hamiltonian from nuclear physics J| a_P JM j and the 
Dicke model from quantum òptics 0, B 13 ■ 

The high degree of similarity between the behaviour 
of these systems suggests an underlying universality, and 
in this paper we explore this universality in terms of a 
quantum mechanical catastrophe theory. 

In its elementary, classical form, catastrophe theory is 
the study of the critical points of potentials, with em- 
phasis on a qualitative understanding of the properties 
of the system as critical points are born, move about, 
merge, and disappear as control parameters are varied 
|lClj |. The best known catastrophe is the cusp, which de- 
scribes the bifurcation of a critical point. The relation 
between entanglement properties and the cusp has been 
noted previously for the Dicke model (TJ, and the impor- 
tance of bifurcations in the appearance of entanglement 
màxima has been conjectured as a general rule Il2|. 
In this paper, we shall explore and expand upon these 
ideas. 

Some properties of the quantum cusp have been dis- 
cussed by Gilmore et al. [iJl but our focus here is differ- 
ent, and the way in which we obtain a quantum model 
from the classical catastrophe differs accordingly. The 
method employed here admits the concept of a macro- 
scopic or semi-classical limit; thus establishing the con- 
nexion with models of quantum phase transitions (QPT). 
The quantum cusp model we construct may be thought 
of as a minimal model that exhibits the salient entangle- 
ment features observed in these models. We study not 
only the cusp, but two further catastrophes — the butter- 
fly and a two-dimensional example — the entanglement 



properties of which expand upon the types of behaviour 
one might expect in more realistic models. 



I. QUANTUM CATASTROPHE MODELS 

We begin by constructing the quantum catastrophe 
models and first consider those derived from catastro- 
phes occurring in a single variable, such as the cusp. 

We take as our model a system of two interacting 
bosonic modes. Let [xi, Pxi) and {x2, PX2) be the (ab- 
stract) position and momentum coordinates representing 
these modes. We assume an interaction between these 
modes such that the interacting system is separable in a 
description in terms of two collective bosonic excitations, 
the coordinates of which we denote {yi,Pyi) and (2/2 , Py2 ) ■ 
We construct the Hamiltonian of one of these collective 
modes yi so that it undergoes the catastrophe. The ques- 
tion that we shall address is then: given the structure of 
the system in terms of the collective modes y, what is 
the entanglement between the original bare modes x? 

We write the Hamiltonian of the collective mode in 
which the catastrophe occurs as 



2ni 



(1) 



with m and lo the characteristic mass and frequency of 
the mode. The potential C/cat(yi) is taken from elemen- 
tary catastrophe theory, and can be written as a power 
series Ucat{y) = X^^i^nJ/i- We rescale the coordinate 
2/1 J/1 y^h/rriLü, and measure the energy in units huj, 
such that 
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which defines the rescaled catastrophe potential Vcatíï/i)- 
Here, /i = moj /Ti is our explícit "macroscopy" parameter, 
which is meant in the sense that the limit /i — > cx) can 
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be thought of either as the limit in which the system 
size (and hence mass m) becomes macroscopic, or as the 
semi-classical limit ?i — > 0. The limit /i — s- cx3 is analo- 
gous to the thermodynamic limit in the QPT models, and 
therein lies the correspondence between these quantum 
catastrophes and the QPT work cited in the introduction. 

The behaviour of the mode described by the Hi is 
largely governed by the fixed points of the classical catas- 
trophe potential 14at(!/i), and this is especially true in 
the limit /i — > (X). By construction the fixed points of 
14at (2/1)7 which we denote y, are of the order y ^ y/Ji, 
and are thus "macroscopic". Expanding Vcat(yi) in Eq. 
(0) about y and taking the limit /i — > 00 we obtain 



H - -l— i — 



yí + v{y). 



(3) 



yi=y 



This effective Hamiltonian describes small 0(1) fluctua- 
tions about fixed point y. The second derivative deter- 
mines the excitation spectrum around the fixed point, 
and V{y) ~ 0(/i) is the energy of the bottom of the 
harmònic potential well in which the system is localised. 
In general, the potential will have more than one fixed 
point and an independent effective Hamiltonian may be 
derived for each. The way in which contributions from 
different fixed points combine to give the overall ground 
state of the quantum system will be treated for individual 
catastrophes. 

The second coUective mode ï/2 is assumed to be simple 
harmònic, and thus the full Hamiltonian of the catastro- 
phe model is 



íícat(y) 



We relate the coordinates of the two coUective modes y 
to those of the bare modes x via the rotation 



yi = cxi + SX2, 
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CX2, 



(5) 



where c = cos(6'/2) and s — sin(6'/2), and 6 refiects 
the degree of mixing. In terms of the x-representation, 
ífcat(x) is not separable, and this rotation generates an 
interaction between the two bare modes x. We quantise 
the coUective coordinates yi and the bare coordinates Xi 
according to 



2/i 



= 2-1/2(61 +Ò,), x,^2-^/\al + a,,), 



(6) 



with momenta defined canonically. In this second 
quantised notation, the two representations are related 
through a two-mode SU(2) squeezing transformation de- 
scribed by the unitary operator W — exp(— |a|a2 -I- 

To make the connexion with a familiar model: the 
above scheme is very similar to the Dicke model in the 
thermodynamic limit. Here, the two bare modes are the 
photon field and the coUective atòmic coordinate, and 
these are related to the coUective excitations (polaritons) 
by just such a squeezing 



In this paper, we consider two one-dimensional catas- 
trophes — the cuspoids and , commonly referred 
to as the cusp and the butterfly. We shall also consider 
a catastrophe that occurs in two dimensions, V^at (2/1,2/2) 
and is non-separable. In this case, we calculate the entan- 
glement between the modes yi and 1/2 with the catastro- 
phe itself providing the interaction between the modes. 
In selecting which catastrophes to study, we require that 
the spectra of the catastrophe be bounded from below 
for all vàlues of the control parameters at finite fj.. 



II. ENTANGLEMENT ABOUT FIXED POINTS: 
/i ^ 00 LIMIT 

For the one-dimensional catastrophes, the two-mode 
Hamiltonian that determines the excitations about yi in 
the ^ 00 limit is 
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(7) 



with = (PV/dyf\^^_^. The ground state wave function 
of the system is thus the Gaussian 
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exp 1 -y2/i 
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which in the x-representation reads 
9 \ 1/4 



*(x) = ( - 
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exp ■{ -^{cxi 
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To find the entanglement of this wave function, we re- 
quire the reduced density matrix (RDM) of one of the 
bare modes, xi, say. This is obtained through p{xi,x'i) = 
/ dx2'9{xi,X2)'i'*ix[,X2) as 



p{xi,x[) = 



Vei(s^ei+c2) 



exp 



\^~a{xl + x[^) + (3xiX2Y;L0) 



where a and /? are coefficients, only the ratio of which is 
important for the entanglement: 



2a _ (ei + 1)2 + 2ei [cot2(0/2) + tan2(0/2)] 
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(11) 



We shall quantify the entanglement in our two mode sys- 
tem with the von Neumann entropy S. The entropy of 
the density matrix p{xi,Xi) is evaluated by comparison 
with the density matrix of a harmònic oscillator at finite 
temperature. Detalls of this approach have been given 
elsewhere 0, and we just give the result here: 



1 ( n / n 

S = <^ — coth — 

log2\2T \2T 



In 
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2T 



(12) 



where the ratio of frequency to temperature of the ficti- 
tious oscillator is given by Í7/T = arccosh(2Q;//?). For the 
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one-dimensional catastrophes, the entanglement is max- 
imised when the squeezing angle is 9 — 7r/2. For this 
choice, Eq. simplifies to 



2a _ ef + 6ei + 1 



(13) 



This procedure is easily adapted to calculatc the entan- 
glement in the two-dimensional catastrophe. 

We now consider our three example catastrophes in 
turn. 



III. CUSP 

The cusp catastrophe, ^+3 is the most familiar and, 
from the point of view of applications, the most im- 
portant catastrophe. With coefficients chosen for con- 
venience, the scaled cusp potential is 



(14) 



We shall only consider a harmònic perturbation here, and 
reserve until later a discussion of the effects of linear per- 
turbations. We also shall set 9 = tt/2 here to give max- 
imum mixing between the modes. This leaves us with a 
single control parameter A. 

The full two-mode Hamiltonian in terms of the creation 
and annihilation operators of the x-modes is 



iï+3(a) 



A 



- — (ajoi + ala2 + 1) 



A-í 



H — («{«2 + aial + a\al + aia2) 

+ -J—{a\ +ai+al + 02)"*. 



(15) 



It may at first appear unusual that the coefíicient of 
ajoí should depend on the parameter A. However, it 
can be shown that, by individually squeezing the coUec- 
tive modes before applying the two-mode SU(2) transfor- 
mation, this dependence on A can be removed. If both 
modes are squeezed identically, the entanglement proper- 
ties of the system are left invariant, since this squeezing 
then represents a global rescaling of the phase space. For 
simplicity though, we retain the form of Eq. H15|l . 

We now consider the fixed points. For A > 0, only 
one stable fixed point exists and this lies at the origin. 
Taking /i — > cx), we see that the excitation energy about 
this fixed point is ei = ^/A. For A < 0, the origin be- 
comes unstable, and two new stable fixed points appear 
at í/1 — fi\A\. In the fj, —^ 00 limit, these two fixed 
points are degenerate and have the same excitation en- 
ergy Cl = 2-y/|yl|. The shape of the potential is sketched 
as insets in Fig. ^ and shows clearly the change of the 
potential from double to single well structure. Note that 
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FIG. 1: Entanglement properties of the cusp catastrophe. 
The von Neumann entropy S in the macroscopic limit fi ^ 00 
(thick line) shows a divergence at the critical value A = 
where the potential changes from a double- to a single-well 
structure (inset sketches) . Numerical results for finite n show 
a peak near this point. Inset (a) shows the scaling with fi of 
the parameter value A* at which the entanglement maximum 
occurs, and (b) shows the value of the entropy 5* at this 
point. 



the form V+3, Eq. (11), is also used in Landau theory 
of phase transition in statistical mechanics, or in quan- 
tum field theory (çi^-model) . Describing the vanishing of 
the excitation energy as ei ~ A^'^ , and the divergence 
of the "correlation length" as ^ = e^^^^ ~ A'^ , we find 
exponents = 1/4 and z — 2. 

We now consider the entanglement. For A > 0, the 
entropy foUows directly from the approach outlined in 
section ^ For ^ < 0, the situation is complicated 
slightly by the existence of two fixed points. With the 
limit ^ 00 taken in correspondence with the thermo- 
dynamic limit, the ground state of the system would be 
an equal mixture of density matrices localised at the two 
fixed points. We prefer here to use the limit ^ —^ 00 
to calculate an approximate wave function for finite but 
large /i. This is obtained by taking a coherent superposi- 
tion of the two localised wave functions and allows direct 
comparison with the numerical results for finite fi. Since 
the two lobes are orthogonal, the reduced density matrix 
of the total system is equal to the sum of the reduced 
density matrices for the two lobes: pi — l/2(p+ -|- p-). 
This is the same result as is obtained if one takes the 
ground-state to be the incoherent mixture; so the dif- 
ference between these two approaches is unimportant. 
However, this will be seen not to be the case when we 
consider the two-dimensional catastrophe. 

From the general theory of entropy 'l6l| we know that 
for p = XiPi with Ai probabilities, the total entropy 
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S{p) is bounded by 



(16) 



In the current situation, since p+ is orthogonal to p_ , the 
upper bound becomes an equality. Furthermore, since 
S{p+) = S{p-), we have S{pi) = Smix + S{p+) with 
S'mix = 1- The mixing entropy represents the contribu- 
tion from the 'global', i.e. macroscopic, structure of the 
wave function, whereas local structure enters through 
the individual S'(p+) terms. If the parity symmetry 
F+3(j/i) = V+3(— t/i) is broken by an additional linear 
term cx yi in the potential, the degeneracy of the two 
fixed points would be lifted and the contribution from 
the mixing entropy S^ix — 1 would disappear. 

The single-well entropy S{p+) is calculated as in sec- 
tion and we plot the total entropy S{p) in Fig. Q 
The similarity between the behaviour of this simple cusp 
model and the QPT models is apparent. At the critical 
point, the entropy diverges as 



S 1 




(17) 



FIG. 2: The von Neumann entropy of the Butterfly catastro- 
phe with A4 = — 4/\/3 as the potential undergoes a double- 
triple-single well transition, both for /i — > cx3 and finite p. 
The profile of the entanglement is very different to that of 
the cusp as the transition here is induced by a level crossing 
in the spectrum. Inset shows scaling of A2 as a function of p. 



i.e., with the correlation length ^, and we thus see "crit- 
ical entanglement" |1 7] ■ 

Numerically obtained results for finite p are shown 
alongside the p ^ 00 result. The value of A for which 
the peak in the entanglement occurs at finite p, A* , scales 
with /i to a very good approximation as A* = cp^ '^^ with 
a numerically determined constant of c = 4.1. This rela- 
tion is plotted in Fig. We mention that the exponent 
of 0.75 «3/4 has been observed numerically for the en- 
tropy in the Dicke model 0- We also investigated the 
value of the entropy S* at its peak (Fig. ^p) but found 
no convincing scaling relation for finite p. 



Case (ii): A4 — — 4/V3. Here we see new behaviour 
absent in the cusp. The A2 parameter range is divided 
up into three regions by the fixed points, 



A2 < 0; y^± 
< ^2 < 4/3; y^O 



y 

A2 > 4/3; y 



y± 

0. 



1/2 



y± 



(19) 



IV. BUTTERFLY 

The second one-dimensional catastrophe that we study 
is the butterfly, which gives rise to the potential 



V+5Íyi 



2 2 



Ai 4 



I 

6^ 



(18) 



The parameter space is two-dimensional (A2,Ai), and 
rather than give a full account of this space, we simply 
look at two representative vàlues of A4 

Case (i): A4 = 0. For A2 > 0, yi — is the only 
fixed point and this has excitation energy ei = ^/A^. 
For ^2 < 0, y = ±.y/I|A2|^/^ are the two stable fixed 
points, both with ei = ·y/2pl2|. Apart from numerical 
coefficients, the behaviour here is the same as that of the 
cusp. This result generalises to all A+k catastrophes: for 
V+fc with Ai = 0,yi > 2 the excitation energy is ^/A^ for 
A2 > 0, and A/(fc - 3)1^2! for A <0, with behaviour like 
that of the cusp. 



Thus, increasing A2 from below zero upwards, the po- 
tential moves through a sequence of first a double, then 
triple, then single well structures, as shown by the insets 
in Fig. El 

The stability or otherwise of the fixed points is only 
part of the story in determining the /i — > 00 ground state 
of the system. For A2 > 4/3 and A2 < 0, the situation is 
straightforward and the ground state is obtained exactly 
as for the two phases in the cusp. In the central region 
< ^2 < 4/3, however, we have three fixed points, and 
their weight in determining the ground state depcnds on 
the energy V{y) of the bottom of the well at y. In the 
p 00 limit, the system will be completely localised in 
whichever of the fixed points has the lowest base energy, 
or, if the energies are degenerate, we take an equal su- 
perposition to describe the large-/i wave function. For 
A2 > 1, y = is the fixed point with lowest energy, and 
for ^2 < 1 the two fixed points at finite displacements 
y — y± have the lowest energy and are degenerate. Only 
at ^ = 1 are all three points degenerate and we have a 
three-lobed wave function. 

This structure is induced by a level crossing in the p — > 
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oo spectrum, with the energy of the double well crossing 
the energy of the single well at ^ = 1. For finite /x, the 
level-crossing is actually avoided, due to the overlap of all 
three wells. This situation therefore bears some snTiilarity 
to that described in Rcf. 18], where a discontinuous 
entanglement was observed at a level crossing associated 
with a first-order QPT. 

Away from the level crossing, the entanglement is cal- 
culated just as for the cusp. In the region of = 1, we 
need to exercise a little care, because the entanglement 
is discontinuous at ^2 = 1. Exactly at this point, the ex- 
citation energies of the three wells do not disappear, but 
rather take the finite vàlues ei = (1, 2, 2). The entangle- 
ment in the central well (with ei = 1) is zero, S'ò — 0, 
since the wave function is circularly symmetric about the 
origin (e2 = 1 as well) and can thus be written as a prod- 
uct state with respect to all co-ordinate systems. The en- 
tanglement for each of the displaced wells is S± ~ 0.197. 
Thus, by combining the appropriate density matrices, we 
íind that for A2 slightly less than unity, the double-well 
state has S — 1.197. For A2 just slightly bigger than 
unity we have S — 0, due to the product state in the 
single well. Directly at ^2 = 1 we have the three-lobed 
wave function, and S = 2/35'+ -I- 1/2S'_ lg3 « 1.716. 
These results plus the corresponding finite fi data are 
shown in Fig. |21 The approach of the finite /i results to 
the /X — > 00 limit is nicely seen, and in particular to the 
limiting value of S w 1.716 at A2 = 1. 

We stress that the entanglement maximum occurs not 
at the value of A2 at which the fixed point becomes un- 
stable, but rather at the level crossing. Moving through 
the points A2 = and A2 = 4/3, where fixed point sta- 
bility does change, nothing special happens to the en- 
tropy (or any other ground-state property), since these 
fixed points do not contribute to the determination of the 
ground state at these vàlues of A2. 

By examining the finite fi data (Fig. we determine 
that the value of A2 at which the entanglement peak oc- 
curs scales as yl* — 1 ~ CQfj,^'^^ with numerical param- 
eters (co,ci) determined to be (—3.55,1.90) to within a 
few percent. 



V. TWO-DIMENSIONAL CATASTROPHE 

The most familiar two-dimensional catastrophes are 
the umbillics with the germs 2/12/2 i yf- However, these 
are unsuitable for our purpose as their spectra are not 
bounded from below and this, in fact, is true of all the 
two-dimensional, elementary catastrophes of Thom |19|| . 
Therefore, we consider the non-simple catastrophe 

Kn = ^Aiyl + vi) + ^{yí + '^ivlvl + vi). (20) 

where we have only included harmònic perturbations as 
before. This catastrophe is described as non-simple be- 
cause the germ (that part proportional to /x~^ in the 




FIG. 3: The von Neumann entropy of the two-dimensional 
molar catastrophe with A = — 1 as a function of 7. Plots of 
the potential for 7 < 1 and 7 > 1 are shown at the top of the 
figure. The origin of the potential is unstable and there are 
four stable potential wells satellite to this. Lower right inset 
shows scaling of 7* as a function of /i. 



above) depends irreducibly on a modulus, 7, whereas 
simple germs have no free parameters. 

The fixed point structure of divides the behaviour 
into three regimes in the fi 00 limit. For A > 0, 
we obtain a single fixed point at the origin, and since 
the ground-state of the system is a product state of two 
Gaussians with the same width, there is no entanglement. 
For A < 0, the origin is unstable; for 7 7^ 1, the system 
possesses four fixed points, as is readily observed from 
the molar-shaped potcntials plotted as insets of Fig. O 
For all 7 > 1, the four stable fixed points lie on the 
lines yi — and í/2 — 0, whereas for 7 < 1 they lie 
on the diagonals yi — ±2/2- In the foUowing, we set 
^2 = — 1 throughout, as the entanglement properties are 
the same for all A2 < 0. We calculate the entanglement 
between modes 2/1 and 2/2 induced by the interaction in 
the catastrophe itself, and do not apply the two-mode 
squeezing. 

We first study 7 > 1 as this is the simpler of the two 
cases. The stable fixed points are given by 

(2/1,2/2) = (±V7^, 0); (2/1,2/2) - (0, íVm)- (21) 

At each fixed point, 2/1 and 2/2 are the excitation coordi- 
nates with excitation energies 

4 = 2; 6^=7-1. (22) 

Excitations in the direction of the displacement are 
described e+. 

The individual wave functions localised around any of 
these fixed points are unentangied, since they are just 
products of Gaussians is the 2/1 and y2 directions. How- 
ever, combining these four functions into the four-lobed 
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wave function that describes the large n limit, the total 
system is entangled. This is solely due to the mixing 
entropy of its four lobed structure. 

We can not calculate the entanglement of this structure 
in the way we did for the one-dimensional catastrophes, 
because the four reduced density matrices of each lobe 
are not orthogonal. This means that the upper bound in 
Eq. H16|) remains as an upper bound, and is not equality. 
Nevertheless, we can proceed as follows. Writing 1 2/1,^2) 
for the wave function of the system localised at (2/1,2/2), 
the four-lobed large-/i wave function can be written as 

I*) = ^{|y,0) + |-y,0) + |0,y) + |0,-y)} (23) 

with y = ^JJi. Given that the individual lobes contribute 
nothing to the entanglement by themselves, we ignore 
their individual structure in this description. In the limit 
/i — > cx), the three single-mode states |0), | ± y) are all 
orthogonal, and thus the RDM of one of the modes pi = 
Tr2|í')(4'| is 

Pi = ^ {( |y) + I - y>) ( (y| + + 2|0)(0|} . (24) 

Furthermore, the orthogonality of these states means 
that this density matrix can be simply treated as a three- 
by-three matrix and the entropy is simply S" = 1, inde- 
pendent of 7 for 7 > 1. 

It is interesting to note that had we taken as the 
ground-state density matrix the incoherent mixture of 
the four contributions, 

P = ^{|y,0)(y,0| + |-y,0)(-y,0| 

-MO,y)(0,y| + |0,-y)(0,-y|}, (25) 

leading to the RDM 

pi = J{|y>(y| + |-^>(-y| + 2|0)(0|} (26) 

and a value of the von Neumann entropy of S = 3/2, 
which is clearly at variance with the numerical results. 

We now consider the region 7 < 1, and for simplicity 
we also assume 7 > 0. The four fixed points are 

where the two ± signs are independent. Each fixed point 
has the excitation energies 

4 = 2; ei = 2^-^. (28) 
1 + 7 

The eigenmodes of the system are not j/i and 2/2, but 
rather lie along, and perpendicular to, the diagonals of 
the í/i-?/2 plane. Each individual fixed-point wave func- 
tion is thus entangled with respect to modes yi and í/2- 
This entanglement can be calculated as in section lïïl 
but here with two excitation energies and the rotation 



between the eigenmodes and the y coordinates. The en- 
tanglement determining parameter 2a /P is evaluated to 
be 

2a _ 4-37^+4^1^ 

f3 ^ 72 ' ^ ' 

from which the single-lobe entanglement follows directly. 

The contribution of the four-lobed structure of the 
large-/x superposition can be assessed as follows. From 
a macroscopic point of view, we can ignore the structure 
of the individual lobes, and write the wave function as 

l*> = l{\y,y) + \y-,y) + \-y,y) + \-y,~y)} 
= {\y) + \^y))^{\y) + \-y))- m 

The second forms clearly shows this wave function to be 
a product state from the macroscopic viewpoint. Thus 
the mixing entropy of forming the four-lobed structure is 
zero, and the entropy of the system is just the single lobe 
entropy above. 

In Fig. Elwe plot these results alongside the numerical 
data for finite /x. The scaling of 7* with /i is observed to 
be 7* — 1 = co/i^^^ with coefficients fitted as (co,ci) — 
(4.93 X 10'', 4.09). 

VI. CONCLUSIONS 

We have constructed and studied a family of quantum 
catastrophe models, and investigated their ground-state 
entanglement properties. The cusp catastrophe, with its 
bifurcating fixed point, demonstrates behaviour that is 
remarkable similar to the QPT models, such as the Dicke 
model — underlining the importance of bifurcations of 
classical fixed points in this context. It should be noted 
that whilst this bifurcation occurs for all vàlues of /i, 
a peak in the entanglement is only observed when /i is 
sufhciently large (/i > 10 here). This illustrates that 
the bifurcation is not, in itself, a sufíicient condition for 
the occurrence of the entanglement maximum, but that 
the system must also be capable of sufficient delocali- 
sation. The butterfly catastrophe displays very different 
behaviour to the cusp — namely a discontinuous entropy 
induced by a level crossing in the macroscopic limit. 

The cusp and the two-dimensional catastrophe demon- 
strate that a mixing term in the entropy can contribute 
to the total entanglement in cases where a wave func- 
tion is split up into localisation areas that are separated 
within (abstract) position space. In particular the two- 
dimensional catastrophe suggests a distinction between 
'global' and 'local' (within the lobes) entanglement, and 
one could speculate that in more complex situations, with 
wave functions split up further and further, a hierarchy 
of entanglement entropies might emerge. 

Our results also have a bearing on the issue of quantum 
chaos and entanglement in such systems, as the model 
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here is capable of emulating the behaviour of more so- 
phisticated nonlinear Hamiltonians, despite being sepa- 
rable — and thus integrable. It is clear that there is no 
unequivocal relation between delocalization and the on- 
set of quantum cl·iaos on one hand and the peaking of 



entanglement on the other. 
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